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Abstract 

Using a Markov chain Monte Carlo approach, we find the allowed parameter 
space of a MSSM model with seven free parameters. In this model universality 
conditions at the GUT scale are imposed on the gaugino sector. We require in 
particular that the relic density of dark matter saturates the value extracted from 
cosmological measurements assuming a standard cosmological scenario. We char- 
acterize the parameter space of the model that satisfies experimental constraints 
and illustrate the complementarity of the LHC searches, B-physics observables and 
direct dark matter searches for further probing the parameter space of the model. 
We also explore the different decay chains expected for the coloured particles that 
would be produced at LHC. 

1 Introduction 

Among the large number of theoretical models proposed to either solve the hierarchy 
problem and/or explain dark matter with a new stable particle, the minimal supersym- 
metric model (MSSM) remains one of the favourite. Supersymmetry not only provides a 
solution to both these problems but also predicts new physics around the TeV scale. The 
main drawback of the MSSM apart from the lack of evidence for supersymmetric parti- 
cles is the large number of unknown parameters most of which describe the symmetry 
breaking sector. With the improved sensitivities of dark matter searches in astroparti- 
cle experiments [H El [HI HJ [5], Ej , the precise determination of the DM relic density from 
cosmology [3 El E] , the latest results from the Tevatron [ini [H] and the precision mea- 
surements, large regions of the parameter space of the supersymmetric models are being 
probed. This will continue in the near future with a number of direct and indirect detec- 
tion experiments improving their sensitivities fTIi [TH| [TH [T^ and most importantly with 
the LHC starting to take data. The LHC running at the full design energy of 14TeV offers 
good prospects for producing coloured supersymmetric particles lighter than 2-3 TeV, for 
discovering one or more Higgs scalars [Ml [I7j and for measuring the rare processes in the 
flavour sector, in particular in B-physics [18]. Furthermore some properties of the sparti- 
cles, in particular mass differences can be measured precisely in some scenarios [IHl E] . 

The first studies that extracted constraints on supersymmetric models worked in gen- 
eral within the context of the MSSM embedded in a GUT scale model such as the 
CMSSM [ini Eni El]. After specifying the fundamental model parameters at the high 
scale, the renormalisation group equations are used to obtain the weak scale particle 
spectrum. This approach provides a convenient framework for phenomenological analyses 
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as the number of free parameters is reduced drastically compared to the general MSSM 
(from O(IOO) to 4 and 1/2 parameters in the case of the CMSSM). The drawback is that 
one is often confined to very specific scenarios, for example in the CMSSM the LSP is 
dominantly bino over most of the parameter space. This has important consequences for 
the dark matter relic abundance. Furthermore it was customary to choose some specific 
values for some of the MSSM or even the SM parameters for a convenient representation 
of the parameter space in two-dimensions. While the link between specific observables 
and allowed region of parameter space is easier to grasp in this framework, the allowed 
parameter space appeared much more restrictive than if all free parameters were allowed 
to vary. 

In the last few years efficient methods for exploring multi-dimensional parameter space 
have been used in particle physics and more specifically for determining the allowed pa- 
rameter space of the CMSSM. This approach showed that the often narrow strips in 
parameter space obtained when varying only two parameters at a time fattened to large 
areas [22l [231 [2H [25] after letting all parameters of the CMSSM and the SM vary in the 
full range. With this efficient parameter space sampling method it becomes possible to 
relax some theoretical assumptions and consider the full parameter space of the MSSM. 
Because the number of experimental constraints on TeV scale physics is still rather limited 
it seems a bit premature to go to the full fledge (9(100) parameters of the MSSM or even 
to the 19 parameters that characterize the model when assuming no flavour structure 
and equality of all soft parameters for the flrst and second generations of sfermions (for 
an approach along these lines see flUi ITTj ). Furthermore many parameters, for example 
those of the flrst and second generations of squarks, once chosen to be equal to avoid 
strong flavour- changing neutral current constraints, do not play an important role in the 
observables selected to flt the model. Here we consider a model where input parameters of 
the MSSM are deflned at the weak scale and we add some simplifying assumptions: com- 
mon slepton masses {Mj^ = Mj^ = Mf) and common squark masses (M^^ = M^^ = Mg 
at the weak scale for all three generations and universality of gaugino parameters at the 
GUT scale. This implies the following relation between the gaugino masses at the weak 
scale, M2 = 2Mi = M3/3. We furthermore assume that At is the only non-zero trilinear 
coupling. While, as we just argued, the flrst assumption should not impact much our 
analysis, the second should certainly be considered as a theoretical bias. This assump- 
tion is however well motivated in the context of models deflned at the GUT scale. Most 
importantly in our approach we keep the higgsino parameter /x and the gaugino mass M2 
as completely independent parameter. The relation between the gaugino and higgsino 
parameters is what determines the nature of the LSP and plays an important role in 
determining the LSP-LSP annihilation in the early universe. In that sense our model 
has many similarities with the non-universal Higgs model which also has fi and M2 as 
independent parameters [28] . 

The observables selected to constrain the model include the relic density of dark mat- 
ter, Qh"^, direct searches for Higgs and new particles at colliders, searches for rare processes 
such as the muon anomalous magnetic moment as well as various B-physics observables. 
Note that the dark matter relic abundance is computed within the standard cosmological 
scenario. The direct detection of dark matter while providing stringent constraint on the 
model introduces additional unknown parameters both from astrophysics and from strong 
interactions. We therefore prefer to consider the direct detection rate as an observable 
to be predicted rather than as a constraint keeping in mind that folding in the astro- 
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physical and hadronic uncertainty could however easily introduce an order of magnitude 
uncertainty in that prediction [29l [30]. 

We find that each individual parameter of the MSSM model is only weakly constrained, 
in particular the parameters of the sfermion sector. The very large allowed parameter 
space only refiects the still poor sampling of the total parameter space by experiments. 
The neutralino sector is better constrained with a preferred value for the LSP of a few 
hundred GeV's and a small likelihood for masses above QOOGeV, similarly charginos above 
1.2TeV are disfavoured. We also find a lower limit on the pseudo-scalar mass as well as on 
tan (3. Furthermore some correlations between parameters of the model are observed, most 
notably the one between /i and the gaugino mass. This is because those two parameters 
determine the higgsino content of the LSP. 

After having determined the allowed parameter space, we examined the predictions 
for direct detection as well as for LHC searches both in the Higgs and SUSY sector as well 
as for B-physics observables. Although each type of search can only probe a fraction of 
the total parameter space we find a good complementarity between the different searches 
with less than 10% of scenarios leading to no signal. For example large signals for direct 
detection are expected in the mixed bino/Higgsino LSP scenario that are hard to probe 
at the LHC. The LHC searches in the SUSY and Higgs sector are also complementary 
and B-observables are specially useful in scenarios with large tan/3 and a pseudoscalar 
that is not too heavy. The predictions for SUSY searches can be different from that 
expected in the constrained CMSSM with in particular a large fraction of models that 
only have a gluino accessible at LHC, the squarks being too heavy. To ascertain how 
experiments that will take place in the near future could further constrain the parameter 
space of the model we consider specific case studies. For example we consider the impact 
of a signal in B[Bs — at Tevatron or of the observation of a signal in direct 
detection experiments. Finally we examine in more details the SUSY signals at the LHC, 
analysing the preferred decay chains for models that have either a gluino or a squark 
within the reach of the LHC. In this analysis we did not include the constraints from 
indirect detection experiments because the rates predicted feature a strong dependence on 
additional quantities such as the dark matter profile or the boost factor. The predictions 
for the rates for p, e"*", 7 will be presented in a separate publication pTj . 

The paper is organised as follows. The model and the impact of various constraints 
are described in section 2. The method used for the fit is described in section 3. The 
results of the global fits are presented in section 4 together with the impact of a selected 
number of future measurements. The SUSY signatures are detailed in section 5. The 
conclusion contains a summary of our results. 

2 Model and constraints 

We consider the MSSM with input parameters defined at the weak scale. We assume 
minimal fiavour violation, equality of the soft masses between sfermion generations and 
unification of the gaugino mass at the GUT scale. The latter leads to M2 = 2Mi = M3/3 
at the weak scale (relaxing this assumption is kept for a further study). We allow for only 
one non-zero trilinear coupling At. For the b-squark the mixing which is oc Af, — /itan/3 
is driven in general by /itan/3 rather than by the trilinear coupling, this approximation 
is however not very good in the small sample of models with /i 100 GeV. Note also 
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that the Higgs mass at high tan j3 can show some dependence on the sbottom mixing. For 
first and second generations of squarks the mixing which depends on fermions masses is 
neghgible except for the neutrahno-nucleon cross section since the dominant contributions 
to the scalar cross section are also dependent on fermion masses. However since the squark 
exchange diagram is usually subdominant as compared to Higgs exchange, the neglected 
contribution of the trilinear coupling falls within the theoretical uncertainties introduced 
by the hadronic matrix elements [30]. Similarly neglecting the the muon trilinear mixing, 
A^, could affect the prediction for (5a^ but this effect is not large compared with the 
uncertainties on the value extracted from measurements. The top quark mass mt is also 
used as an input although it has a much weaker influence on the results than in the case 
of GUT scale models. For the latter the top quark mass enters the renormalization group 
evolution and can have a large impact on the supersymmetric spectrum in some regions 
of the parameter space. In the general MSSM the top quark mass mainly influences the 
light Higgs mass. We fix as{Mz) = 0.1172 and mb = 4.23 GeV. The free parameters of 
our MSSM model with unified gaugino masses, MSSM-UG, are 

fi,M2,Mi,Mg,At,ta.np,MA,mt (1) 

The range examined for each of these parameters is listed in Table [H MSSM-UG has a far 
more restricted set of paramters than the general MSSM, still this model will show how 
the possibilities for SUSY scenarios open up. The observables that will be used in the fit 
are listed in Table [21 We first review the expectations for the role of each observable in 
constraining the MSSM parameter space. 



Table 1: Range of the free MSSM-UG parameters. 



Symbol 


stands for 


General range 




parameter 


[-3000, 3000] GeV 


Ms 


Gaugino mass, 2Mi = M2 = M3/3 


[30, 2000] GeV 


Ml 


Common slepton mass, Mj = ^ 


[50,4000] GeV 




Common squark mass, Mg = rriq^ ^ 


[50,4000] GeV 


At 


Trilinear coupling of t 


[-3000, 3000] GeV 


tan/9 


tan [3 


[5,65] 


Ma 


Mass of CP-odd Higgs boson 


[100,2000] GeV 




mass of t- quark 


[165, 180] GeV 



The most powerful constraint is VLh"^ . Since in the MSSM with gaugino universality 
there are four different mechanisms for efficient neutralino annihilation each calling for 
a different combination of SUSY parameters, each individual parameter is weakly con- 
strained when exploring the full parameter space of the model. The main mechanisms for 
neutralino annihilation are: 

• annihilation of a bino LSP into fermion pairs, this necessitates light right-handed 
sfermions. 
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• annihilation of a mixed bino/higgsino into W pairs, this means /i ~ Mi. Some 
contribution from chargino and neutralino coannihilations is also expected. 

• annihilation near a (light or heavy) Higgs resonance, a LSP with a non-zero higgsino 
component is required but the higgsino fraction can be very small if m^o mh/2. 
For heavy Higgs exchange this process is more efficient at large tan/? due to the 
enhanced couplings of the heavy Higgs to b-quarks and r leptons. 

• Coannihilation with sfermions, this usually means those of the third generation since 
they are lighter due to mixing. Recall that we are assuming only one slepton mass 
Ml and one squark mass Mq. 

The lower limit on the Higgs mass from LEP strongly constrains the small tan/? 
region. The upper range for the top quark mass and/or a large mixing in the stop sector 
are favoured. The latter means large values of \At\. 

The muon anomalous magnetic moment still indicates a Sa deviation from the SM 
expectations [32] • The SUSY loop contributions can lead to 5a^ > when sparticles are 
light. This is the case for the chargino/sneutrino loop and the neutralino/smuon loop 
when n > 0. The latter is enhanced at large values of tan/?, in this case it is easier 
to accomodate 6a^ > even with a SUSY spectrum at TeV scale. Recall that we have 
neglected A^, although the mixing is usually dominated by /xtan/?, this introduces some 
restriction on the parameter space. 

The branching ratio for B{b — > s'-f) depends on the squark and gaugino/Higgsino sector 
as well as on the charged Higgs. A light pseudoscalar mass is permitted only if squarks are 
light as well. The squarks then partially cancel the large pseudoscalar contribution but 
only for yU > 0. Furthermore when /i > 0, the chargino contribution is negative relative 
to the SM one so B{b 57) can drop too low for light charginos. This is particularly 
true when At is large and negative, meaning a large mixing in the stop sector. 

The process B{Bs f^~^f^~) receives SUSY loops contributions from Higgs exchange 
as well as from chargino and sfermions. In particular the amplitude for Higgs exchange 
is enhanced as tan^ /?. The largest contribution are therefore expected for large tan (3 and 
light Ma, particularly for light charginos and sleptons. 

In the MSSM, the branching ratio for B^ — > is suppressed relative to the standard 
model prediction by the charged Higgs contribution. This contribution is enhanced as 
tan^ /?. We will use the observable R{Bu rur) which gives the ratio between the SUSY 
and SM predictions. Due to large uncertainties in hadronic matrix elements the standard 
model branching fraction is not known with a good precision, so this observable is not as 
powerful as other B-physics observables to constrain the parameter space of the MSSM. 

LEP limits on the light Higgs and on sparticles constrain the chargino and neutralino 
sector as well as the charged sfermions. We have not included the recent Tevatron limits 
on the squark and gluino mass as well as the limits on the chargino mass from the trilepton 
search. When presenting our results we will comment on the impact of the new Tevatron 
results. 

^Note that the trilepton search in pp — s- Xi^X2 extends the LEP constraint on the chargino mass only 
when sfermions are light, no constraint is found in the large scalar mass limit (large mo) at least in the 
context of the CMSSM. The large fraction of our models which have heavy squarks will therefore not be 
probed by the Tevatron even if they feature light neutralinos and charginos. 
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If (chain==l) 



chain=0, frb4-+ 




fra=frb 





^ If (count < Nmax) ^ 




Initialize Xa, Xb, La, Lb 
fra, frb, Oa. Ob, count,chain 



Calculate Spectrum 





count++ 


Get Xb with proposal 
function Xb=Q(X,Xa) 







Figure 1: Details of the MCMC algorithm used in this paper. Xa and Xb are the points 
in the parameter space, Oa and Ob are the values of the observables at points Xa and Xb 
respectively, La and Lb are corresponding likelihoods and fra and frb are the frequencies 
at two points. 

3 Parameter sampling method 

In this section we attempt to motivate the MCMC method in a heuristic way, for a detailed 
treatment see Ref. [23] 

The likelihood L is the probability distribution function (PDF) p{d\m) for a data 
set d being reproduced by an assumed model m. In our case, we assume the model 
m =MSSM-UG and the data set to be reproduced is given in Table O This is a top- 
down approach where by varying the model parameters we find a suitable m to maximise 
p{d\m). However, in a bottom-up approach one would like to know the probability of 
model m being correct once given the data d, i.e. p{m\d). From Bayes' theorem we have 

»(m) 



p{m\d) = p{d\m) 



p{d) 
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where p(m) is the (absolute) probabihty of the model m being correct and p{d) is the 
total probability of reproducing the data d for all possible models. The p{d) is hard to 
estimate since it requires knowledge of all possible models. Thus, the absolute value of 
p(m|(i) is hard to estimate. However, it is possible to estimate the relative correctedness 
of a model mi against another model m2 given the data d by taking the ratio 

p{mi\d) p{d\mi) p{mi) 
p{m2\d) p{d\m2) p{m2)' 

For simplicity, we take p{mi) = p{m2), that is the probability of two models being correct 
is same or we chose a fiat prior over the model space. In this paper this translates to 
using a flat prior over the parameter space of MSSM-UG as listed in Table [B Thus with 
flat prior we have 

p{mi\d) L{mi) 
p{m2\d) L{m2)' 

This clearly hints that if we sample the parameter space using a directed random walk 
with transition probability proportional to min(l, L(mi)/L(m2)), the sampling density 
will be proportional to the likelihood ratio or the relative correctedness of the model. The 
details can be found in [33] and a representation of the algorithm used is given in Fig. [H 
There are two points in order: For the sampling density to be proportional to the 
ratio of the hkelihoods the proposal PDF, (^(X, Xa), should be symmetric about Xa, 
where X, Xa are point in the parameter space. We use the Gaussian distribution in each 
parameter direction and the proposal PDF looks like: 



2 (AX 



^2 



The second point is about the definition of likelihood. The likelihood L is a product of 
likelihood for each of the observables. For the observable with definite measurement the 
likelhood function is a Gaussian 



G(0,Oe.p,AO)=exp 



2 (A0)2 



Here, Oexp is the central value of the observable, AO is the la error and O is the value 
of the observable at the proposed point X in the parameter space. For observables with 
upper/lower limits the likelihood function is a smooth step-like fuction 

F(0,Oe.p,AO) = - ^ 



l + exp[±{0-Oe.p)/AO] 

The +ve sign is for the lower limit and the — ve is for the upper limit. Here Oexp is the 
95% exclusion limit and AO is about 1% of Oexp to roughly emulate the 95% exclusion 
limit. The list of observables used to calculate the likelihood is given in Table [2] alongwith 
the correponding experimental values or limits. 

4 Results 

We have first scanned linearly over the 8 parameters of the MSSM-UG using the Markov 
Chain Monte Carlo (MCMC) method just described. The particle spectrum was computed 
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by SoftSusy2.0 [M] and fed to micrOMEGAs2.2 [Ml ESI ED] for the computation of all 
DM observables as well as of constraints on the parameters of the supersymmetric model. 
The SM value for the branching B{b ^ 57) calculated in [35] was shifted in order to 
match the NNLO result of [37]. 

The range chosen for all parameters is listed in Table [T] and the six observables used in 
the fit listed in Table [21 In addition we have imposed the LEP limits on sparticles , 
mi...) as defined in micrOMEGAs and for the likelihood we have assumed a binary step 
function. We then examined the allowed regions for each input parameter as well as the 
predictions for physical parameters and observables. We have not imposed the Tevatron 
limit on gluino and squarks > 308 GeV and > 379 GeV [38]. However after fitting 
the model we have checked a posteriori that these constraints did not affect much our 
analysis. Because of the universality condition the gluino mass limit is always satisfied 
while the squark limit is not satisfied in less than one per-mil of our scenarios. 



Table 2: Observables used in the fit. 



Observable 


Limit 


Likelihood function 




(27.8 ±8.5) X 10-1° 


G(a:,27.8 x 10-i°,8.5 x 10-^°) 


B{h S7) 


(3.55 ±0.24) X 10"^ 


G(x,3.55 X 10-^0.24 x 10"^) 




0.113 ±0.0105 


G(x, 0.113, 0.0105) 


B{Bs ^ 


< 0.8 X 10-^ 


F(x,0.8 X 10-^-0.8 X 10-9) 


R{B T+v) 


1.11 ±0.52 


G(x, 1.11,0.52) 


rrih 


> 114.5 GeV 


F(a;, 114.5,0.6) 




171.4 ±2.1 GeV 


G(x, 171.4,2.1) 



Table 3: Range of the MSSM parameters. 



Parameter 


68% C.L. 


95% C.L. 


/X 


197.7 


1193 


-983.5 


2471 




259.8 


1077 


119.6 


1845 


Ml 


398.4 


2270 


225.4 


3700 




1824 


3602 


1236 


3938 


At 


-1735 


2239 


-2785 


2876 


tan (3 


19.5 


65.0 


8.15 


65.0 


Ma 


537.1 


1489 


370.2 


1904 




169.6 


173.7 


167.6 


175.6 
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Figure 2: Likelihood function for the model parameters, from top left: 
fi,M2,Mj,Mq, tan Ma, all dimensionfuU parameters in GeV. 



4.1 The allowed parameter space of the MSSM 

We found wide allowed regions for each model parameter, the distributions are displayed 
in fig. [5] and the 68% and 95% limits in Table [31 

• As expected, > is preferred, this is mainly due to the B{h 57) and 5a ^ 
constraints. Note however that one can get a reliable global fit even with 5a^ = 
so that the 95%C.L. extends to the negative /i region, see Table [H 

• The gaugino mass reaches almost the maximum value probed, M2 < 1845 GeV at 
95%C.L. Charginos above the TeV scale are somewhat disfavoured (see the 68% C. 
L. at 1.1 TeV) because they give little contribution to 6a^, furthermore the light 
neutralinos LSP annihilate more efficiently except when they have a large Higgsino 
component or can annihilate near a resonance. 

• The slepton mass distribution is peaked at 500GeV with a long tail that extends 
to almost the upper limit of the region scanned. The sleptons just above the LEP 
limit are disfavoured with 225GeV < M; < 3700 GeV at 95%C.L. This shows that 
although light sleptons can contribute to the annihilation of a bino LSP and to 5a^, 
present data can be accomodated with heavy sleptons. 

• Heavy squarks are preferred with Mq > 1.24 TeV at 97.5%C.L. Quarks below the 
TeV scale tend to give too large corrections to B-physics observables and in partic- 
ular to B{b 57). Although the quark contribution can be compensated by that 
of a light pseudoscalar, this requires fine-tuning and has a small likelihood. There 



9 




Figure 3: 2d-likelihood functions for a) Mi — fi b)M2 — /i c)Ma — M2 showing the 
68%C.L.(dark grey) and 95% C.L. (hght grey). All parameters in GeV. 

is no upper limit on the squark mass as it can reach the upper limit of the range 
used in the scan, 4 TeV. 

• Light pseudoscalars are disfavoured with Ma > 370 GeV at 97.5% C.L. While it 
is possible to have a good fit to the data with lighter pseudoscalar masses this 
occur only after fine tuning the sfermion masses in order to satisfy both 6a^ and 
B{b 57). Therefore these models have a small likelihood. 

• The distribution for tan/5 is skewed towards the upper range and tan/3 > 8.15 at 
97.5% C.L. Large values of tan/3 make it easier to escape the LEP bound on the 
Higgs mass, facilitate annihilation of the LSP through Higgs exchange and can also 
help explain 6a ^. 

• There are no preferred value for At although the distribution is skewed towards 
At > 0, see Table [3l 

A few correlations among input parameters are observed. These are driven to a large 
extent by the relic density constraint. First fi and Mj are correlated: a large /i requires 
light sleptons while heavy sleptons need fi not much above the TeV scale. When sleptons 
are heavy the LSP must have a higgsino component to annihilate efficiently. Conversely 
when /i is large and the LSP is bino, light sleptons are needed for efficient annihilation into 
fermion pairs. Second M2 is also correlated with /i and M4: when M2 exceeds 1.2TeV, 
one needs either ^ M2/2 ^ Mi or M2 — Ma- The first correlation is strong and 
means that a heavy LSP with a significant higgsino component annihilates efficiently. 
The second correlation corresponds to m^o ^ Ma/^ with the LSP annihilating near a 
heavy Higgs resonance. This correlation is weak as we have just mentionned, other large 
values of M2 are allowed when the LSP has a large higgsino component. 

4.2 The particle spectrum 

From this allowed parameter space we can determine the favoured region for the masses 
of new particles that can be searched for at LHC. 

• The light Higgs is SM like and rrih = 120.2+^9 GeV at 95%C.L. 
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Table 4: The 68% and 95% allowed range for sparticle masses and for the higgsino fraction. 



Mass 


68% C.L. 


95% C.L. 


TD I, 




1 99 6 


115 3 


1 94 8 


771 TT 




1 489 


370 


1 900 

X C/ WW 


Tfl -0 
Xi 




^99 


^9 9 


874 


771 -n 


1 OA 




1 04 


1 1 


771 -n 

A3 






1 




771 -n 


41 9 






9^1 n 


771 -L 




74Q 


1 n"^ 


1 1 


771 -L 

^xt 




1 (S4n 


9(^9 


94Qn 


Tfl ~ 


41 n 


99^0 


zou 


^^71 n 


7T1 ~ 


406 


9970 




3700 

O 1 WW 


TD ~ 


1 03 


9960 


1 03 


3690 

OW C/W 


Til ~ 


457 


9980 


976 


371 

O 1 xw 


771 ~ 


1 R4n 


OUtJU 


1 9'^n 


Oc/c/ W 


771 ~ 


1 R^O 


3660 


1 940 


4000 

TiW WW 


7T1 7 


1 ^^n 


0\J ( u 


1 940 


401 

^W 1 w 


771 ~ 


1 ^^0 




1 940 


4000 

t:W WW 




1 ^90 




1 9nn 


Oc/UW 


7T1~ 


1 8^0 




1 950 


C/W 


Til 7 
"Hi 


1770 


3590 


1130 


3930 




1880 


3660 


1300 


4000 


m-g 


952 


3320 


478 


5400 


fn 


1.58 10-3 


0.297 


3.39 10-^ 


0.659 



• The LSP mass is in the range 52 GeV < m^o < 873 GeV at 95%C.L. The lower 
bound results from the assumption of GUT scale gaugino mass universality and the 
LEP limit on charginos. The peak in the distribution around 55 GeV corresponds 
to LSP annihilation near a light Higgs resonance. 

• The lightest chargino (as well as X2) could be just above the LEP reach and lies 
below m-+ < 1182 GeV at 97.5%C.L. There is a peak in the distribution near 
lOOGeV which is correlated with the LSP peak correponding to annihilation near 
the light Higgs. 

• The gluino mass is related to the neutralino mass via the universality condition, the 
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97.5%C.L. lower limit gives rrig > 477 GeV and the maximum can exceed 5 TeV. 
However we have Mg < 3.3 TeV with 84% C.L., which is slightly beyond the reach 
of the LHC with maximum luminosity. 

• As mentionned before, squark masses are in the TeV range. Because of the univer- 
sality assumption on squark masses the mixing in the stop sector implies that the 
stop is the lightest squark. We find m^^ > 1130 GeV. Squark masses can reach all 
the way to the upper end of the range of the scan, 4TeV. There is no guarantee that 
squarks are within the reach of the LHC. 

• Slepton masses reach almost all the way to 4TeV. Since masses of sleptons and gaug- 
inos are uncorrelated, a large fraction of scenarios have mj^^ > m^o,m-+ (around 
75% of allowed scenarios). For these scenarios the preferred decay channel of X2 
is three body and has a branching fraction around 3% into each flavour of charged 
lepton. In this model all slepton masses are identical except for mixing effect so 
only one mass is displayed in Fig. H] 

The higgsino fraction of the LSP strongly influences the properties of the DM, fn = 
l^isl"^ + |A'^i4p where the neutralino mixing matrix is defined in the basis Xi = NuB + 
N12W + N13H1 + Ni4^H2. The higgsino fraction spans a wide range 3.4 x 10~^ < fu < 0.66 
at 95%C.L. The distribution is peaked around 0.15 < fn < 0.35. A small higgsino fraction 
is found either for light LSP's and for LSP's annihilating near a Higgs resonance. The 
higgsino fraction has an impact on the spectrum and in particular on the mass difference 
m-+ — m^o. A large mass difference is expected for a small higgsino fraction (since in that 
case m^+/m^o f» 2) whereas when fn > 0.15, m-+ — m^o < Mw and m^o — m^o < Mz- 
This means that for a large number of scenarios that satisfy this condition the decay 
channels xf — > Xi^^ X2 ~^ Xi^ are forbidden and only 3-body decays are possible since 
sleptons are often heavier than the chargino. 



Table 5: Range of predictions for observables 



Observable 


68% 


C.L. 


95% 


C.L 






0.102 


0.123 


0.091 




0.133 


Sa^ X 10^° 


1.49 


27.2 


-1.45 




37.5 


B{B X7) X 10^ 


3.19 


3.62 


3.08 




3.86 


B{Bs X 10^ 


1.82 


3.31 


1.48 




20.4 


R{B Tv) 


0.698 


0.981 


0.382 




0.997 


4i (pb) 


2.51 X 10-1° 


6.76 X 10-« 


9.77 X 10-12 


2.24 X 10-^ 



4.3 Other observables 

Form the allowed parameter space we then examine the predictions for the observables 
used in the fit together with the direct detection cross-section. 
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Figure 4: Distributions for sparticle masses: m^o , m -+ , m^o , m -+ , rrig, m/i,me^,m^^ and 
for the higgsino fraction, f^- There is a scale factor of 10^ on the y-axis. 

• The observables Qh'^ and B{h ^ 57) are well distributed around the central values 
used in the fit. As expected, B{h ^ 57) lies near the upper end of the allowed range 
for small Ma-, /i, M2 and and are enhanced at very large values of tan (3. 

• The branching ratio B{Bs is peaked around 2.5 x 10~^ although a long 
tail extends to 2. x 10~^ at 95%C.L. Values near the upper limit are found for 
Ma < 600 GeV, tan/5 > 50, for not too heavy charginos, m-+ < 750 GeV and for 
a LSP with a small higgsino fraction, fn < 0.05. 

• At 97.5%C.L., R{B tu) > 0.38, as for other B-physics observables, the largest 
suppression to R{B tu) are found for small MA,fi,M2 and Mj as well as for 
tan/3 > 50. 

• The muon anomalous magnetic moment extends over more than an order of mag- 
nitude and in particular can include negative values. The peak of the distribution 
lies near 10~^°, much below the central experimental value we have used indicating 
that our model does not provide a good fit to this observable. Large deviations in 
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5a ^ are found for light sleptons and neutralinos/charginos, this is why a better fit to 
the data is found for low values of M2 and M^. This is ilUustrated in Fig. [7K,b. This 
observable also plays a role in setting the upper limit on the neutralino/chargino 
masses. 




Figure 5: Distributions for B{h sj), B{Bs — > ), (5a^, cr^p. The y-axis is in units of 
10^. 



Note that if Sa^ turns out to be consistent with the SM prediction, as indicated by 
some preliminary results [M], the distributions for M2, Mj would shift to larger values. The 
dependence of the B-observables on some of the most relevant parameters are displayed 
in Fig. [7b-i, in particular the enhancement with the light pseudoscalar mass for B{Bs — > 
fi~^fi~), B{b — > s'j),R{B — > Tu). In fact a heavy pseudoscalar leaves little prospects for 
observing large deviations in B-observables, Ma > 1 TeV implies B{Bs f^^f^~) < 
5 X 10"^ and R{B tu) > 0.70. Similarly small values of tan/3, for example tan/5 < 20 
implies B{Bs — > < 2.5 x 10~^ and R{B ru) > 0.85. The potential enhancement 

at large values of tan/3 is illustrated only for B{Bs — > fi^fi~). The largest values for 
B{Bs — > /i are also found at low values of M2, they are therefore correlated with 
large deviations in 6a^ [ID]. 
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Figure 6: Predictions for elastic scattering cross sections as a function of rn^o, /i, Ma and 
fn- Default values of micrOMEGAs for the quark content of the nucleon are taken. On 
the top-left plot the CDMS-II exclusion is shown by a solid line, while the dashed line 
show the factor 10 improved limits. 



The predictions for the neutralino nucleon elastic scattering, a^^, vary over 5 orders of 
magnitude. The distribution peaks at values near the experimental upper limit, see Fig. [5] 
and can even exceed the present limit. The LSP mass dependent CDMS limit is displayed 
in Fig. [6^ [2] . The light Higgs exchange in general dominate SI interactions. Therefore 
large cross sections are expected when the LSP has a large coupling to the light Higgs, 
this means some higgsino component. Fig. ^jp,d. For example a LSP with a higgsino 
component fn > 0.2 necessarily implies a^p > 3.2 x 10~^ pb. These values will be probed 
in the near future with for example XenonlOO [12]. On the other hand a pure bino LSP has 
a cross section at least two orders of magnitude below the experimental limit. The second 
Higgs scalar as well as squarks can also contribute significantly when they have a mass 
comparable to the LSP. There is however no direct correlation between the direct detection 
rate and Ma- Note that the enhancement of the heavy Higgs contribution relative to the 
light Higgs expected at large values of tan 13 is tamed because the pseudoscalar is always 
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much heavier than the hght Higgs (recall that > 370 GeV). 




Figure 7: 2D contours for some observables, from top left to bottom right: a) Sa^ — Mj 
b) Sa^ - m^o c) B{Bs - d) 5(5, - At e)B{Bs - tan/? 

f) BiBs ^ - Ma g)B{B, - h)B{b ^ 37) - Ma i) i?(5 ^ rz/) - M^. 

The 68%C.L. is in dark grey and 95%C.L. in light grey. 
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4.4 Complementarity of different LHC searches and direct dark 
matter searches 

From the above discussion, it is clear that B-physics observables, (5a^, direct detection 
as well as direct particle searches at LHC are sensitive to different sectors of the MSSM. 
Direct searches at LHC are especially powerful to probe the coloured sector while heavy 
Higgs searches probe the large tan j3. Direct detection probes best the higgsino LSP. 

To illustrate quantitatively the complementarity of the various new physics signals at 
LHC and in direct detection we compute the fraction of models in the allowed parameter 
space that lead to a signal in 4 different channels: direct searches of sparticles at LHC, 
searches for heavy Higgs at LHC, deviation in B{Bs and DM direct detec- 

tion. Here we choose B{Bs — > /x'''/^") as a representative observable from the B-sector 
to simplify the discussion. The same exercice with another B-observables would lead to 
a similar conclusion. The signal in each channel is define as follows. We assume that 
coloured sparticles are within reach at LHC-14TeV if their mass is below 2 TeV. This cor- 
responds roughly to the reach for gluinos when squarks are heavy and C = 60 fb~^ [T6lfT7]. 
This value for the squarks reach is conservative but we use it for simplicity, recall that 
with a high luminosity, £ = 300/6~\ coloured sparticles can be discovered up to 3 TeV. 
To define the discovery reach for the heavy Higgs we use the results of the CMS study 
which gives the discovery region in the tan/5 — Ma plane for C = 30/6"^ ^7\. This is 
based on the associated production of a Higgs with b-quarks, gg, qq —>■ bbh with the Higgs 
decaying into tau pairs. In the MSSM, this process is accessible only at large values of 
tan j3 because of the tan /3-enhanced couplings of the heavy Higgs to bb, tt pairs. For 
example the pseudoscalar mass can be probed up to Ma < 800 GeV for tan/5=50. 

For the process Bg we assume the value for the branching B{Bs fJ'~^fJ'~) < 

5x 10"^ which will be probed even with a low luminosity at the LHC. Finally for the elastic 
scattering cross section our definition of a signal corresponds to one order of magnitude 
improvement over the current CDMS limit p]. With various experiments increasing their 
sensitivities regularly, this level should be reached in the near future [12]. 

Table 6: Prospects for Higgs searches, direct detection and B{Bs f^~^fJ'~) in scenarios 
with heavy coloured sparticles, rriq, nig > 2 TeV 



Heavy Higgs at LHC 


pb) 


B{Bs ^ /i+/i") 
> 5. X 10-9 < 5. X 10-9 


Yes 
24.2% 


20.6% (> 10-9) 
3.6%(< 10-9) 


5.6% 15% 
1.7% 1.9% 


No 
75.8% 


54.9%(> 10-9) 
20.9%(< 10-9) 


0.3% 54.6% 
0.7% 20.2% 



First we consider the case where sparticles are too heavy to be accessible at LHC, 
in these scenarios it is clearly essential to consider alternative signals of new physics 
and/or dark matter. In fact a large number of the allowed scenarios (almost 29%) predict 
a heavy coloured spectrum. We will refer to this as set A in the next section. This 



17 



is easily explained by the fact that no observable (except weakly for 5a^) requires a 
super symmetric contribution as long as the DM can give the proper relic abundance. The 
best alternative at LHC for a sign of physics beyond the standard model would be the 
search for a heavy Higgs d with nearly 24.2% of scenarios predicting a signal while few 
deviations from the SM in B{Bs — > fi~^fi~) are expected (only in 5.6% of the total number 
of scenarios in set A). Furthermore signals in B-physics are mostly found in scenarios 
for a light pseudoscalar Higgs and very large values of tan (3 where one also has a heavy 
Higgs signal. The complementarity between the direct detection and collider searches is 
clearly evident in Table El Nearly half (54.9%) the models lead to a signal only in direct 
detection experiments and about 20.6% of models predict a signal in both the Higgs and 
the direct detection. Conversely if no signal in direct detection is observed in the near 
future the model would be severely constrained (with 24.5% of scenarios remaining). Of 
course this statement should be modulated by the fact that there are large theoretical 
and astrophysical uncertainties in the prediction of the neutralino proton scattering cross 
section. Combining all the channels we are left with a small fraction of models with no 
signal in SUSY or Higgs searches at colliders, B physics or direct detection, the fraction 
corresponds to 20.2% of the subset with heavy sparticles which means about 6% of the 
sample of allowed scenarios. We have also checked that there are no other SUSY signals at 
LHC in these scenarios with heavy squarks and gluinos. We have computed the trilepton 
signal for direct neutralino/chargino production, pp — ^ XiXi — > 3/ + and found 

a maximum cross section of a few fb at ^/s = 14 TeV. The reasons for the small cross 
section are the heavy squarks and the heavy neutralinos and charginos, the (97.5%C.L.) 
lower limits are m^o.m^o.m^o > 285,366,371 > GeV and m-+ > 354 GeV. The direct 

Al' X2' A3 ' ' Xi 

production of sleptons is also below the LHC reach for SOfb"^ since the lower limit on 
slepton masses are , mg^ > 300 GeV. 

Table 7: Prospects for Higgs searches, direct detection and B{Bs f^^fJ'~) in scenarios 
with light coloured sparticles, rriq < 2 TeV and rrig < 2 TeV 



Heavy Higgs at LHC 


pb) 


B{Bs ^ /i+/i-) 
> 5. X 10-9 < 5. X 10-9 


Yes 
24.5% 


21.9% (> 10-9) 
2.7%(< 10-9) 


7.2% 14.7% 
1.6% 1.1% 


No 
75.5% 


53.7%(> 10-9) 
21.8%(< 10-9) 


0.6% 53.1% 
1.1% 20.7% 



The results for the case where sparticles are accessible at LHC leads to roughly the 
same conclusion, see Table [71 One difference is a small increase in the fraction of models 
that predict an observable deviation in B{Bs — > fi~^fi~). This is particularly true if squarks 
below the 2TeV scale are present. In the case where a heavy Higgs and squarks are 
accessible at LHC, more than 35% of the models also predict B{Bs IJ'^ > 5. x 10-9. 

^Here we have not included the possibiUty of discovering SUSY particles in Higgs decays 
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4.5 Impact of improved sensitivity 

Rather than just looking at the discovery potential we also examine the impact of a 
signal in the near future in direct detection or in B{Bs fi^fi~). First let us comment 
on the impact of the present CDMS limit on direct detection. Recall that because of the 
astrophysical and nuclear uncertainties we have not used this observable in the fit. We 
applied this constraint a posteriori on the selected models and looked at the impact on 
the SUSY spectrum. The most noticeable effect of this constraint is that it removes some 
models with 65 < m^o < 170 GeV especially when those models have a large value for 
the higgsino fraction fn- Models with m^o ^ 50 GeV that have a small higgsino fraction 
are not affected. Furthermore some of the models with low values of and large tan f3 
are disfavoured. To have an indication of the impact of a signal in direct detection we 
then considered the case of a signal a^^ = 1. x 10"*^ pb and allowed a factor 3 theoretical 
uncertainty. With such a measurement one could constrain the allowed parameter space 
although no specific correlation with other observables are observed, only a lower bound 
on the LSP higgsino fraction is found. 

We have also considered the impact of a signal at the Tevatron in B{Bs fi'^fi'), say 
B{Bs — s> fi~^fi~) > 1.8 X 10~® which represents the ultimate reach of Tevatron [41]. We 
found that in this case a lower bound on tan/5 > 32 with 97.5%C.L could be extracted 
(though a tail extends till tan/3 ~ 18) and furthermore that a heavy Higgs signal at 
LHC would almost be guaranteed (in 94.6% of the scenarios). Finally the models with 
a signal within the Tevatron reach in B{Bs f^~^fJ'~) all predict a small cross section in 
the trilepton channel {(jJif^ < 0.1 fb) even before applying cuts. This points towards a 
complementarity between the two channels as was found in the MSSM [1^, although as 
we have mentionned before only a small fraction of the MSSM-UG models are within the 
Tevatron reach for the trileptons. 

5 LHC forecasts 

At last we examine in more details the SUSY signatures at the LHC. For this we split the 
scenarios allowed at 95% C.L. in different sets according to which sparticles are within 
the reach of the LHC. We have already mentionned that a significant fraction (28.7%) of 
our scenarios predict a coloured spectrum that is just too heavy (set A). Another large 
fraction (30.5% ) corresponds to the case where the gluino is the only coloured sparticle 
that could be discovered at LHC (set B) while a squark is also within reach in 22.0% (set 
C). In 19% of our scenarios (set D) only squarks are accessible. The main difference from 
the CNMSSM is that in the latter the squark masses receive an important contribution 
from M3 due to the renormalization group running so the cases where only a gluino is 
accessible at LHC are confined to the focus point region. The precise fraction of models 
that have sparticles within reach of the LHC do depend on the range on the parameters 
used in the scan. In particular it is clear that a wider range for the squark masses would 
have lead to a larger fraction of models with squarks too heavy to be produced at the 
LHC. 

For two sets of scenarios we then examine the main branching fractions of sparticles 
concentrating on squarks, gluinos and on the neutralino/charginos produced in their decay. 
Our goal is to point out the different decay channels available. We do not attempt to 
analyse the feasibility of extracting signals over background in our scenarios, this is beyond 
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the scope of this paper. 

Note that in some scenarios, the chargino/neutrahno can be produced directly leading 
to a signal in the trilepton channel, pp X^xt — Sl+E^^^"^. We do not here analyse direct 
production specifically but we have checked that the Tevatron only weakly constrained 
our allowed scenario, only about 1 per-mil of our models have a trilepton cross section 
a{3l + E™"^) > 0.0.08 pb before applying cuts. 

5.1 Scenarios where rrig < 2TeV, Mq > 2 TeV. 

In the MSSM-UG the scenarios where only gluinos are accessible at LHC represent a large 
fraction of the scenarios that provide a good fit to the data. To examine the characteristic 
decay chains in these scenarios, we have, for each point in our MCMC chain, computed 
the branching ratios of gluinos into neutralino/charginos as well as the 2- and 3-body 
branching ratios of X2 X3 with micrOMEGAs. The decay chains involving a chargino 
can be dominant but we concentrate on the heavy neutralinos because the leptonic decay 
mode provide a clear signature and the kinematic endpoint in the lepton invariant mass 
distribution further allow for a determination of the mass difference of the heavy neutralino 
with the LSP. 

In these scenarios the squarks are heavy so the gluino decays only via 3-body, g 
xtff or g ^ XaJI' leading to final states with many jets as well as missing energy 
from the LSP. The branching fraction into third generation quarks is particularly relevant 
since it has been shown HI] that requiring tagged b-quark jets in the final state 
helped reduce the SM background and thus could extend the LHC reach in this channel. 
In our sample scenarios, the total branching fraction for 3-body decays of gluinos into 
third generation quarks {t and b) B{g xQQ') = ^gQQ varies over a wide range, from 
0.1 — 0.8 and features two peaks around 25% and 70%, see Fig. [HkEI The X3 channels 
each contribute below 20%, see Fig. [8b. The larger branching fractions Bgqq are found for 
low values of see Fig. [9l this is mainly because in this case all charginos and neutralinos 
can be produced in gluino decays and contribute significantly (especially charginos). The 
branching fraction BgQQ also increases with tan this is because of the enhanced coupling 
of the higgsino to bb. 

Next we have computed the decay modes of the X2 Xs- The X2 decays dominantly 
via 3-body, the branching fraction of neutralinos into each fermion pair is typical of 
the Z decay (around 3% for each leptonic mode), see Fig. [TUl In some scenarios the 
two-body leptonic channels are accessible, in this case the X2 either decay exclusively 
into the ligthest slepton (the stau) or into all sleptons including the invisible mode into 
sneutrinos, see Fig. [TOb . Finally the decay channels into a gauge or Higgs boson and the 
LSP dominate for sufficient mass splitting between the neutralinos, the dominant mode 
is X2 ~^ Xih unless only the Z is kinematically accessible. For Xg a larger fraction of 
scenarios have 2-body decay modes, in particular the decays into neutral and charged 
gauge bosons, Xs ~^ Xi^^oy Xi'^- Decays into sleptons are as before confined to the 
scenarios with light sleptons. Fig. [TTl 

Finally we have compiled the fraction of models (among set B only) that feature 
each dominant decay mode, the results are presented in Table [HI The largest sample 
(35%) corresponds to both neutralinos decaying into 3-body channels. In these scenarios 

■^For reference recall that in a typical focus point scenario the branching ratios of gluinos into third 
generation quarks reach 72% [15] . 
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Figure 8: Distribution of the gluino branching fractions in third generation quarks ei)BgQQ 
b) Br{g — > x^QQ). There is a scale factor of 10^ on the y-axis. 
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Figure 10: Distribution of X2 branching fractions for the channels a) Xi^~^^ {f^~^t^ ) b) // 
c) Xi^; Xih- There is a scale factor of 10^ on the y-axis. 

jj, < 350 GeV, /i — Ml is small and the LSP has a significant higgsino component so that 
the direct detection rate is large (a^p > 7 x 10~^ pb). The scenarios where Xs instead 
has 2-body decays also have in general a LSP with an important higgsino component 
and a large a^^ except when Xs ~^ Xi^~ is the dominant channel. Then we found 
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Figure 11: a) Distribution of Xs branching fractions for the channels a) Xi'C'e b) 
// c) XiZ, Xih,XiW- There is a scale factor of 10^ on the y-axis. 



Table 8: Dominant decay mode of X2 Xs in models with rrig < 2 TeV and Mg > 2 TeV 
and fraction of allowed models corresponding to each dominant decay mode. 





xl 


Fraction of models 


xVf 


xUf 


35.2% 




xfw^ 


4.2% 


xVf 


XiZ, Xih 


8.6% 




11 


0.2% 


x'lZ 


xtw^ 


6.0% 


Xlh 


X^Z (or xfW^ ) 


14.4% 


11 


11 


4.6% 




xfw^ 


6.1% 


xlff 




3.8% 


xlz 




2.8% 


x> 




6.7% 


11 




5.0% 



predominantly m^o ^ 60 GeV with a bino LSP annihilating efficiently through a light 
Higgs resonance. 

A large number of scenarios (31% of set B) correspond to both X2 X3 having 
2-body decay modes. The case where the dominant decay mode is X2 ~^ XiZ correspond 
to a small value for /i — Mi, this means a dominantly higgsino LSP. In general though 
the dominant modes are X2 ~^ Xi^ and Xs ~^ XiZ unless sleptons are light. In nearly 
19% of sceanrios B, is large and the Xs which is higgsino cannot be produced in gluino 
decays. The X2 decays predominantly via 2-body channels, //, xih or xiZ or into 3- 
body X2 Xlff with the leptonic mode offering a robust signature. Typically in these 
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scenarios the elastic scattering cross sections is small apart from the cases where the heavy 
Higgs is light enough to contribute significantly. 

In conclusion, although the focus point like scenario with gluinos decay with large 
fraction into b quarks and leptonic 3 body decays of X25 X3 constitute the largest sample, 
the channels with two body decays of neutralinos are also very important. Among these 
cases where leptonic decays of X2 dominate constitute 16% of our sample while those 
where the only final states include Z(W) or h constitute about 30% of our sample. Note 
that the gluino only scenario could be quite challenging for the LHC especially if the 
gluino is near 2TeV. We have checked that the scenarios that would be accessible with a 
lower luminosity, say nig < 900 GeV for C = 30fb~^, share the main features we have just 
described. Obviously because of the lighter spectrum a smaller fraction of scenarios have 
X3 accessible in gluinos decays. Furthermore the fraction of scenarios with 2-body decays 
is smaller with for example hardly any scenario where both X2 X3 decay via 2-body. 

5.2 Scenarios where Mq < 2 TeV, rrig > 2 TeV 

In the MSSM-UG a small fraction of allowed scenarios (set D) predicts that only squarks 
are lighter than 2TeV. In this case, the main contribution to SUSY production is direct 
production of squark pairs. In the MSSM-UG, the masses of sleptons and squarks are 
not correlated ,so sleptons are not always light enough to be produced in neutralino 
decays. The conventional SUSY searches which involve a decay chain with a neutralino 
Q ~^ ~* ~^ Xi^^j will be superseded by the two-body decays, X2 ~^ x\Z, Xi^ or 
xfW or if kinematically forbidden by 3-body decays. 

We have computed the branching fractions for squarks, taking as an example the ul 
and ur. We have then looked for the dominant mode into either a chargino, a heavy 
neutralino or the direct decay into the LSP. In each case we have also examined the decay 
mode of the neutralinos and charginos that occur in the squark decays. The higgsino 
fraction of the LSP is an important factor that determines the dominant decay mode. 
We therefore analyse the dominant decay modes separately for the case of a bino, mixed 
or higgsino LSP. Note that a mixed LSP constitutes more than half of our sample of 
scenarios. 

The right-handed squark decays mainly into the bino component which means that 
the dominant decay is in general ur —>■ uxi, leading to only jets and missing energy. The 
decay ur —>■ mx-j can occur only when the LSP is a higgsino {fn ^ 50%) so that the 
Xi channel is suppressed. In this case however the heavier neutralino decay mostly via 
3 body decays into the lightest chargino or neutralino -|- jets leading again to signatures 
with jets and missing energy. 

The left-handed quark which couples strongly to the wino and/or higgsino component 
has a wide variety of decay modes. The frequency of each dominant decay chain are 
displayed in Table [9] for each LSP configuration. For the bino LSP the dominant mode is 
usually Ul ^ dxf with typical branching fractions around 60%. The chargino will decay 
either into XiW or //' when light sleptons are present. The subdominant mode in those 
scenario is ul ux^ with xH ~^ Xi^? X^h- The decay chains are similar to those of 
the CMSSM. In some cases the second chargino, a mixed higgsino/ wino, is kinematically 
accessible and the dominant mode will be — > dxt with subdominant decays into x\'^ 
and Xid- xt decay preferentially into X2^ 1^ other neutralinos as well as into II. 
The Higgs can be produced in either xt ~^ xt^ further in X2 Xi^- A fraction of 
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Table 9: Dominant decay chains for ul for a bino, mixed or higsino LSP. 





bino 
fn < 0.01 


mixed 
0.01 < fH< 0.5 


higgsino 
fH > 0.5 


Ut — > U 

A 1 


7.7% 


2.5% 




A3 Al 




2.0% 


0.3% 
0.3% 
7.4% 


UL xtd xt ^ X\W 
xt 11' 

Al 


12.4% 
8.9% 






Ut — > Y^C? — * y'?!^ 

"^ij A2 A2 A? ^ * 

xl ^ xth 

xt ^ 


2.7% 
1.1% 


37.2% 
1.0% 
12.5% 


2.5% 
0.6% 
0.6% 


Fraction of sample 


32.8% 


55.3% 


11.9% 



models (7.7%) feature the dominant decay into the LSP ul wx?- Because the squark 
ql has a suppressed rate to the bino, this channel is dominant only when other two-body 
channels are kinematically forbidden. 

For a mixed LSP (0.01 < /h < 0.5) the relative importance of the various decay 
channels shifts. The decay ul ^ uxi is dominant in less than 3% of the cases although 
because of the higgsino component of the LSP this can occur even when heavier neutralinos 
are kinematically accesssible. By far the most frequent dominant decay is ul — > dxt with 
significant branching fractions in ux^ or dxt- The heavier chargino always has two-body 
decay modes, xt ~^ Xa^ (preferably X2^) or xt ~^ xt^- The X25 X3 and xt turn 
feature mostly 3-body decays. Note that decay modes into Higgs bosons xt ^ xt^ can 
involve even the heavy Higgs bosons. As usual when light sleptons are present the decay 
xt — > II' can be dominant. 

For the higgsino LSP (/// > 0.5), the dominant mode is either x^u or xtd with some 
contributions from xtd and x^u. The xt channel has similar decay chains as the mixed 
LSP except that the dominant mode is usually xt xti^ rather than channels involving 
Higgses. The Xs can in a few cases decay via two-body, xt"^ or //, but in most cases it 
decays via three-body dominantly into X3 xtff- These decays mainly give signatures 
into jets and missing energy. The xt produced in squark or neutralino decays will also 
decay via three body final states. 

In summary qi decays dominantly into heavy charginos with further decay chains 
involving other chargino/neutralino states. Decay chains involving slepton production 
dominater only in 25% of scenarios. Finally recall that the elastic scattering cross section 
also differs significantly depending on the nature of the LSP giving an opportunity to 
correlate SUSY signals at LHC with those of direct detection. For the bino LSP, cr^p < 
10~^ pb while a^p > 10~^(5 x 10~^) pb for the mixed (higgsino) LSP. 
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We do not discuss in detail the case where both squarks and gluinos are below 2TeV. 
The decay chains can be rather complicated with the possibility of producing the gluino 
in squark decay and vice-versa. The case where the gluino is heavier than the squarks 
features the same decay chains for the squarks as the case just discussed. 

6 Conclusion 

Increasing the number of free parameters as compared to the CMSSM model has opened 
up the possibilities for supersymmetric scenarios that are compatible with all experimental 
constraints and this even maintaining the universality of gaugino mass. Although the 
parameter space of the model is still not very well constrained, we found that the most 
favoured models have a LSP of a few hundred GeV with a significant higgsino fraction 
(> 10%). Contrary to the CMSSM case the higgsino LSP is not fully correlated with 
a very heavy squark sector although all our scenarios favour squarks above the TeV 
scale. A very light pseudoscalar is also disfavoured with Ma > 370 GeV, this means 
that large deviations from the SM in B-physics observables are expected only in a small 
fraction of allowed scenarios. Our favoured scenarios predict few signals at the Tevatron, 
the pseudoscalar Higgs as well as the coloured sector are too heavy to be accessed by 
direct searches. Only very few scenarios have a potentially large enough rate for trilepton 
searches at the Tevatron. The complementarity between future experiments to probe 
this class of models was emphasized. Even though SUSY or heavy Higgs signals are not 
guaranteed at LHC, the majority of allowed models predict at least one signal either at the 
LHC (including the flavour sector) or in future direct detection experiment. Furthermore 
the light Higgs is expected to be around 120GeV with SM-like couplings. 

We have also explored the various dominant decay chains for gluinos and squarks that 
could be produced at LHC in the MSSM-UG as well as for the heavy neutralinos appearing 
in the decays of these coloured sparticles. We found that for models with gluinos accessible 
at LHC, a significant fraction of the heavy neutralinos produced decayed dominantly 
into a gauge or Higgs boson. Furthermore states which decayed into sleptons are rarely 
dominant. We also showed how the preferred squarks decay channels are determined to a 
large extent by the neutralino composition. Whether one can exploit these decay chains 
to determine some properties of the sparticles remains to be seen. In our analysis the 
relic density measurement plays the dominant role in constraining the model, since the 
relic density computation implicitly assumes a standard cosmological scenario, relaxing 
this requirement affects significantly the allowed parameter space of the model. 

Finally we comment on the difference between our results and other recent analyses 
done within the framework of the MSSM with 24 parameters, either using a MCMC 
likelihood approach or applying 2a constraints on each of the observables [26l HE] . First 
these studies were done in a more general model than the one we have considered, with in 
particular no universality condition on the gaugino masses. This means that the LSP can 
have a significant wino component and therefore is more likely to be at the TeV scale as 
was found in [27] using linear priors. Recall that a TeV scale wino annihilates efficiently 
into gauge bosons pairs. The analysis of [27] also emphasizes the prior dependence with 
a generally much lighter spectrum using log priors. This is due mostly to the poorly 
constrained parameter space [TTl HH] . As in our analysis squarks and sleptons ran over 
the full range allowed in the scan and the pseudoscalar mass can be very heavy. 
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The analysis of used a different statistical treatment but most importantly did not 
require that the neutalino explained all the DM in the universe (only an upper bound on 
Qh"^ was imposed). This means that a large number of models with small mass splitting 
between the LSP and the NLSP appeared in the scan calling for a careful study of collider 
limits. In our approach such models are ruled out since they have Qh^ ^ 0.1. This analysis 
further emphasized the light SUSY spectrum in their scans so naturally found preferred 
LSP mass below the TeV scale. 
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